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Abstract — We propose a parametric model of the saturated 
Permanent-Magnet Synchronous Motor (PMSM) together with 
an estimation method of the magnetic parameters. The model 
is based on an energy function which simply encompasses the 
saturation effects. Injection of fast-varying pulsating voltages 
and measurements of the resulting current ripples then permit 
to identify the magnetic parameters by linear least squares. 
Experimental results on a surface-mounted PMSM and an 
interoir magnet PMSM illustrate the relevance of the approach. 

Index Terms: Permanent magnet synchronous motor, mag- 
netic circuit modeling, magnetic saturation, energy-based mod- 
eling, cross-magnetization 

I. Introduction 

Sensorless control of Permanent-Magnet Synchronous Mo- 
tors (PMSM) at low velocity remains a challenging task. 
Most of the existing control algorithms rely on the motor 
saliency, both geometric and saturation-induced, for extracting 
the rotor position from the current measurements through 
high-frequency signal injection lEl- However some mag- 
netic saturation effects such as cross-coupling and permanent 
magnet demagnetization can introduce large errors on the 
rotor position estimation Q, ||4l. These errors decrease the 
performance of the controller. In some cases they may cancel 
the rotor total saliency and lead to instability. It is thus 
important to correctly model the magnetic saturation effects, 
which is usually done through d-q magnetizing curves (flux 
versus current). These curves are usually found either by finite 
element analysis FEA or experimentally by integration of 
the voltage equation (5], (6\. This provides a good way to 
characterize the saturation effects and can be used to improve 
the sensorless control of the PMSM El, O. However the FEA 
or the integration of the voltage equation methods are not so 
easy to implement and do not provide an explicit model of the 
saturated PMSM. 

In this paper a simple parametric model of the saturated 
PMSM is introduced (section [III); it is based on an energy func- 
tion ||9l, ifTOll which simply encompasses the saturation and 
cross-magnetization effects. In section Hill a simple estimation 



method of the magnetic parameters is proposed and rigorously 
justified: fast-varying pulsating voltages are impressed to the 
motor with rotor locked; they create current ripples from 
which the magnetic parameters are estimated by linear least 
squares. In section |IVl experimental results on two kinds of 
motors (with surface-mounted and interior magnets) illustrate 
the relevance of the approach. 

II. An energy-based model for the saturated 
PMSM 

A. Energy-based model 

The electrical subsystem of a two-axis PMSM expressed in 
the synchronous d — q frame reads 



dt 



dt 



^^d + -J79q 

dt 
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(2) 



where 



are the direct-axis flux linkages due to the 



current excitation and to the permanent magnet, and (j)q is the 
quadrature-axis flux linkage; Ud^ Uq are the impressed voltages 
and idjiq are the currents; is the rotor (electrical) position 
and R is the stator resistance. The currents can be expressed 
in function of the flux linkages thanks to a suitable energy 
function l-L{(t)dAq) by 

id = din{(j)dAq) (3) 

H = d2n{^dAq). (4) 

where dkl-L denotes the partial derivative w.r.t. the k^^ variable, 
see 13, Col; without loss of generality H(0, 0) = 0. 
For an unsaturated PMSM this energy function reads 

wliere Ld and Lq are the motor self-inductances, and we 
recover the usual linear relations 

id = di'H{(f>d, 



d2n{4>d,<P,) = j^- 



Notice the expression for 1-L should respect the symmetry 
of the PMSM w.r.t the direct axis, i.e. 

-(/),) =H((/)^, (/),), (5) 

which is obviously the case for Hi. Indeed, O-© is left 
unchanged by the transformation 

(^d' ^d^ ^'d^ ^q^ ^q^ '= i^d^ Ud, id, -(t>q, -Uq, -Iq, S)] 

this implies 



I.e. 



Therefore 



d{H{(l)d, ■ 



K) 



diH{(l)dAq) 



dU 

d(l)d 



dU 

d(j)q 



i>d, -(j^q) dll-L{(t)d, -(t)q) 
= din{(j)d,(l)q) 

d(pd 

{(t^d, -(t^q) = -d2U{(t)d, -(t^q) 
= d2H{(j)d,(t)q) 

-rry9d,(Pq)' 



Integrating these relations yields 

H{(j)d, -(l>q) = y^{(l>d, (l>q) + Cd{(t)q) 
1-L{(l)d, -(t>q) = 1-L{(t>d, (t>q) + Cq{(j)d), 

where Cd^ Cq are functions of only one variable. But this makes 
sense only if Cd{(t)q) = Cq{(j)d) = c with c constant. Since 
H{0, 0) = 0, c = 0, which yields ©. 

B. Parametric description of magnetic saturation 

Magnetic saturation can be accounted for by considering 
a more complicated magnetic energy function having Hi 
for quadratic part but including also higher-order terms. From 
experiments saturation effects are well captured by considering 
only third- and fourth-order terms, hence 

l-L{(t)dAq) = y^MdAq) 

3 4 

i=0 i=0 

This is a perturbative model where the higher-order terms ap- 
pear as corrections of the dominant term 1-Li. The 9 coefficients 
aij together with L^, Lq are motor dependent. But (|5]) implies 

<^2,i = <^o,3 = cif3,i = cifi,3 = 0, SO that the energy function 
eventually reads 
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(b) (f)q(id = Constant, iq) 
Fig. 1. Flux-current magnetization curves (IPM) 

From ©-d?]) and © the currents are then explicitly given by 

id = dil-L{(l)d,(l)q) 

d 2 2 3 2 

-l^d 



d2H{(l>dAq) 

-A 

- + 2ai,2^(i0g + 2a2,2^^0(? + 4ao,4^g, 



(8) 



ao,4^ 



q' 



(6) 



which are the flux-current magnetization curves. Fig. [T] shows 
examples of these curves in the more familiar presentation of 
fluxes w.r.t currents obtained by numerically inverting (IS)-©; 
the motor is the IPM of section [iVl 

The model of the saturated PMSM is thus given by dU)-© 
and ©-©. It is in state form with as state variables. 

The magnetic saturation effects are represented by the 5 
additional parameters as^o, <^i,2, <^4,07 cif2,2, <^o,4- 



C. Model with id-, iq as state variables 

The model of the saturated PMSM is often expressed with 
id-tiq as state variables, e.g. f5|. Starting with flux-current 
magnetization curves in the form 



(t>q = ^qiid^iq) 

and differentiating w.r.t time, O-© then becomes 

i-ddK^d.lq)-^ + ^^dq\}d,^q)-^ = Ud - Rid + ^ 



(9) 
(10) 



Lqdiid, h)-^ + Lqq{id, H)-^=^ 

where 

fLdd{id,iq) Ldq{id,iq)\ ^ fdi^d{id,iq) d2^d{id,iq) 

\Lqd(id,iq) Lqq{id,iq)) \dl^q{id,iq) 

Though not always acknowledged Ldq and Lqd should be 
equal. Indeed, plugging ©-(111) into (l9l)-([TQb gives 

^d = ^d{din{^d.^q).d2n{^d.^q)) 
^q=^q{din{^dAci).d2n{^dAq))- 

Taking the total derivative of both sides of these equations 
w.r.t. (j)d and (})q then yields 

1 0\ ^ (LdddiiU + Ldqdi2n Lddd2in + Ldqd22n\ 

I J \Lqddlin^Lqqdl2n ^,^^21^ + ^,,922^; 

Ldd Ldq\ fdiiU d2iV} 

Lqd LqqJ \di2H 822^^^ 

Since di2l-L = d2iH the second matrix in the last line is 
symmetric, hence the first; in other words Ldq = Lqd. 

To do that with the model of section III-BI the nonlinear 
equations dU)-® must be inverted. Rather than doing that 
exactly, we take advantage of the fact the coefficients ai^j are 
experimentally small. At first order w.r.t. the ai^j we obviously 
have (j)d = Ldid + 0{\ai^j\) and (j)q 
Plugging these expressions into d?])-© we easily find 



Lqiq + 0{\ai^j\). 



t>d = Ld{id - 3a3,oi^^^^ - ai,2Ll 



AaAQL^d^^ 



d^d 



2a2,2LdLqidiq) 



(11) 



Lq{iq - 2ai^2LdLqidiq - 2a2,2LlLqi']iq 



- 4ao,4L^g + r). (12) 

Finally, 

Ldd{id, iq) = Ld{l- Qas^oLdid - 12^4,0^^^^ - 2^2,2^^^^) 

LdqiU, iq) = Lqdiid, iq) = -2LrfL^Zg (^1,2 + 2a2,2^d^(i) 
Lqq{id,iq) = Lq{l - 2^1,2^^^^^ - 2^2,2^^^^ " 12^0,4^^^^) 



III. A PROCEDURE FOR ESTIMATING THE MAGNETIC 
PARAMETERS 

A. Principle 

To estimate the 7 magnetic parameters in the model, we 
propose a procedure which is rather easy to implement and 
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Fig. 2. Experimental illustration of equation dlSl : time response of ^ 



reliable. With the rotor locked in the position ^ = 0, we inject 
fast- varying pulsating voltages 



Ud{t) 
Uq{t) 



Ud + Udfi^t) 
Uq^Uqf{nt), 



(13) 
(14) 



where Ud^Uq^Ud^Uq^Vt are constant and / is a periodic func- 
tion with zero mean. The pulsation ^ is chosen large enough 
w.r.t. the motor electric time constant. It can then be shown, 
see section IIII-C[ that after an initial transient 



id{t)=~id^idF{^t)^0{^) 
i,{t) 



la + 



(15) 

i^F{nt)^0{^), (16) 

where id = ^jiq = ^^idjiq are constant and F is the 
primitive of / with zero mean (F has clearly the same period 
as /); fig. [21 shows for instance the current id obtained for 
the SPM of section [IVl when starting from id{0) = and 
applying a square signal Ud with Q = bOOHz, Ud = 23^ 
and Ud = 30F.J3n^the other hand using the saturation model 
the amplitudes id^iq of the current ripples turn out to be 



Id 



ft\Ld 

+ 12a4,oLlilud 



2a2,2Lqiq{2LdidUq + LqiqUd) 

Gas^oLdidUd + 2ai^2LqiqUq 

\ / ll _ _ _ 

n KIT' ^ '^^'^aLdidi2LqiqUd + LdidUq) 

-f 12ao,4i^g^g2q + 2ai^2{LdidUq + LqiqUd)^. 



(17) 



(18) 



As id-,iq can easily be measured experimentally, these ex- 
pressions provide a means to identify the magnetic param- 
eters from experimental data obtained with various values 

of Ud,Uq,Ud,Uq. 

B. Estimation of the parameters 

Since combinations of the magnetic parameters always 
enter (fTTt-JTSl) linearly, they can be estimated by simple linear 
least squares; moreover by suitably choosing Ud^Uq^Ud^Uq, 
the whole least squares problem for the 7 parameters can be 
split into several problems involving fewer parameters: 



with Ud 



0, hence id 



0, and Ud = 



(resp. Uq = 0) equation (fTTl) (resp. equation ([TS]) ) reads 

with Uq = 0, hence iq = 0, and 2g = 0, (fTTl) reads 

- + Q(^3,oLdid + 12a4,o^d^d 



Ud 



(19) 



(20) 



Notice (fTSl) reads = hence provides no information 
with = 0, hence id := 0, and Uq = 0, ([T7l)-([T8]) read 



1] 



with i^fi 



0, hence id := 0, and i^c^ 

- 2Un ^ -T 



■aioi 



i -^(1- 



12ao,4Lqiq 



(21) 
(22) 

0, (O-dll]) read 
(23) 
(24) 



(resp. Lq) is then immediately determined from ([T9] ): 
aa^o and a4,o are jointly estimated by least squares from (|20] |: 
<^2,2, <^i,2 and ao,4 are separately estimated by least squares 
from respectively (EB, and (l24l) . 



C. Justification of section \III-A\ 

The assertions of section IIII-AI can be rigorously justified 
by a straightforward application of second-order averaging 
of differential equations ||TT1 p. 40]. Indeed the electrical 
subsystem O-© with locked rotor (i.e. ^ = 0) and input 
voltages ([T3l)-([T3l) reads when setting r = Vtt 

^ ^^^"^^ ^df{r) - Rid{(i)d, (t^q)) 



a^q ± 

d7 = ?^^"« 



-Uqf{T) - Riq{(t)d,(t)q)). 



(25) 
(26) 



This system is in the so-called standard form for averaging, 
with a right hand- side periodic in r and ^ as a small 
parameter. Therefore its solution is given by 



where (0^, is the solution of the system 



(27) 
(28) 



^ = ud-md{^l^l) 



dt 



Uq - Riq{(tPd^(t)l) 



obtained by averaging the right-hand side of (1251)- (|26I). After 
an initial transient (^^(r), ^^(r)) asymptotically reaches the 
constant value (^d-, 4^q) determined by Ud = Rid{4>d^ 4^q) and 

Uq = Riq{4)d,4>q)- 



Plugging ([27l)-([28l) with t = ^ into dT])-®, and expanding 
along powers of ^ then yields 

, , - F{nt) fud - ^ - ^ 

id{t) = H 7^ — -r- + 6a3,o0d^d + 2ai^2(PqUq 

\J^d 

+ 12a4,o^i2rf + 2a2a{2^d^qUq + 4>lud) ) +0(^) 



F(nt) fUq ^ , - ^ - ^ , 

■ ( — + 2ai^2[(pdUq + (pqUd) 



+ 2a2,2(20d0giid + (l)^dUq) + 12ao,40^^g I + 0{^) 



where id 
to express 



id{(i>d,(l>q) and 



= iq{4>d^4^q) There remains 
iq in function of idjiq- Rather than exactly 
inverting the nonlinear equations (IT])-®, we take advantage 
of the fact the coefficients aij are experimentally small. At 
first order w.r.t. the aij we have 0^ = Ldid + and 
(j)q = Lqiq + 0{\aij\). Using this in the previous equations 
and neglecting 0{-^) and 0{\aijf) terms we eventually 
find ([T5])-([T8]). Using directly ([TTJ-dni) yields of course the 
same result. 



IV. Experimental Results 

A. Experimental setup 

The methodology of section [Till is tested on an interior 
magnet PMSM (IPM) and a surface-mounted PMSM (SPM) 
with rated parameters listed below. The setup consists of an 
industrial inverter with a 400V DC bus and a AkHz PWM 
switching frequency, 3 dSpace boards (DS1005 PPC Board, 
DS2003 A/D Board, DS4002 Timing and Digital I/O Board) 
and a host PC. The measurements were sampled also at AkHz. 
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Pole pairs 


6 
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Rated power 


200 W 


1200 W 


Rated current 


1.2 A 


3.4 A 


Rated speed 


1800 rpm 


400 rpm 


Rated torque 


1.06 N.m 


29 N.m 


Resistance 


12.15 ^ 


6.69 ^ 



B. Experimental results 

With the rotor locked in the position 6> = 0, a square wave 
voltage with frequency Vt = bOOHz and constant amplitude 
Ud or Uq (30V for the IPM, 40V for the SPM) is appHed 
to the motor. But for the determination of LdjLq where 
Ud = Uq = 0, several runs are performed with various Ud 
(resp. Uq) such that id (resp. iq) ranges from —2A to +2^4 
with a 0.3A increment (IPM), or from —SA to SA with a 
0.5 A increment (SPM). The estimated parameters are listed 
below; the uncertainty in the estimation stems from a ±10mA 
uncertainty in the current measurements. 
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Fig. 3. Measured values (circles) and fitted curve (solid line) for (|20). 





IPM 


SPM 


Ld (mH) 


91.9 ±5 


155.4 ± 10 


Lq {mH) 


45.8 ±1 


58.6 ±2 


a3,o {A.Wb-'') 


7.70 ±0.11 


5.01 ±0.11 




5.35 ±0.61 


4.83 ±0.27 


0^4,0 (A^^"^) 


19.42 ± 1.34 


1.83 ±0.28 


a2,2 (Aiy^-^) 


22.18 ±2.80 


8.76 ± 1.03 


ao,4 (AW6-^) 


6.62 ± 0.42 


1.18 ±0.17 



The good agreement between the fitted curves and the 
measurements is demonstrated for instance for (|2Qt on Fig. [3] 
and for (l22l) on Fig. IH Notice (l20l) illustrates saturation on a 
single axis, while (1221) illustrates cross-saturation. 

C. Validation 

The estimation procedure relies on (|2Q]|-(|24]|, with either 
Id 7^ or , i.e. current vectors with angles 0°,90°, 180°, 270°. 
To check the validity of the model tests were conducted with 
current vectors with various angles and magnitudes on the 
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Fig. 4. Measured values (circles) and fitted curve (solid line) for (22] . 



whole operating = a/^^ ± ranging from OA to 2 A with a 
0.3A increment for the IPM, and from OA to 5.5A with a 0.5 A 
increment for the SPM). Fig. [5] shows for instance the results 
for a 60° current angle; there is a good agreement between 
the measured values and those predicted by the model. 

As a kind of cross-validation we also examined the currents 
time responses to large voltage steps. Fig. [6] shows the good 
agreement between the measurements and the time response 
obtained by simulating the model with the estimated parame- 
ters; it also shows the differences with the simulated response 
when the saturation effects are omitted. Fig. [7] shows the 
good agreement also between the "measured" flux values (i.e. 
obtained by integrating the measured currents and voltages) 
and the flux values obtained by simulation. 

V. Conclusion 

A simple parametric magnetic saturation model for the 
PMSM with a simple identification procedure based on high- 
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Fig. 5. Measured values (circles) compared to model-predicted values (solid 
line) for a 60° current angle. 



frequency voltage injection have been introduced. Experimen- 
tal tests on two kinds of PMSM (IPM and SPM) demonstrate 
the relevance of the approach. This model can be fruitfully 
used to design a sensorless control scheme at low velocity. 
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